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We study the energetics of island formation in Stranski-Krastanow growth within a parameter-free 
approach. It is shown that an optimum island size exists for a given coverage and island density if 
changes in the wetting layer morphology after the 3D transition are properly taken into account. 
Our approach reproduces well the experimental island size dependence on coverage, and indicates 
that the critical layer thickness depends on growth conditions. The present study provides a new 
explanation for the (frequently found) rather narrow size distribution of self-assembled coherent 
islands. 

PACS numbers: 68.65.+g, 68.55.-a, 81.10.Aj 
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The surface morphology of overlayers in heteroepitax- 
ial growth has attracted intense interest because of its 
importance for basic science and applications in optoelec- 
tronic devices. Experiments |l|^ showed that heteroepi- 
taxy in systems with a lattice constant difference > 2%, 
such as InAs/GaAs |-§, Ge/Si (l) and InP/InGaP [|, 
follows the so-called Stranski-Krastanow growth mode 
three-dimensional (3D) dislocation- free (so-called co- 
herent) islands form on top of the wetting layer. These 
small coherent islands are often found to have a very nar- 
row size distribution and are promising to be used 
in quantum dot light emitting diodes (LEDs) and lasers. 

It is commonly agreed that the energetics of strain re- 
lief plays a key role in the growth process: Islands form, 
instead of a uniformly strained, epitaxial film, because 
the gain of elastic relaxation energy in an island over- 
compensates the cost due to the increased surface energy 
by islanding. It is tempting to attribute the observed 
island size distribution to a minimum of the free energy 
of the system. However, an equilibrium theory with only 
two energetic contributions, a positive one from the is- 
land surface energy {E ~ V'^^^, V is the quantum dot 
volume) and a negative one from the elastic relaxation 
energy (E ^ V), fails to predict a finite equilibrium size 
of the islands. Instead, the energy gain from strain re- 
lief always prevails for sufficiently high coverages, ren- 
dering larger islands more stable than smaller ones. In 
order to cope with this difficulty, several additional ef- 
fects, e.g. contributions from intrinsic surface stress or 
from interactions between islands 0,^, have been in- 
voked. Priester and Lannoo proposed a mechanism 
in which 2D platelets act as precursors for the forma- 
tion of 3D coherent islands, thus determining their size. 
Most recently, the observation of island ripening has 
made it doubtful if the islands can be interpreted at all 
as structures in total equilibrium. 

In this Letter, we show that the narrow size distribu- 
tion of the coherent islands can be understood as the 
result of the system being trapped in a constrained equi- 
librium state, where the size is determined by the island 



density and the nominal coverage. In the constrained 
equilibrium theory, the existing nuclei grow to a size de- 
termined by the energetic balance that governs material 
transport between the wetting layer and the islands. This 
allows us to derive an optimum islands size for a fixed 
coverage and island density from a parameter-free ap- 
proach. The elastic energy in both the islands and the 
substrate is calculated within continuum elasticity the- 
ory. Using density-functional theory within the local- 
density approximation, accurate surface energies are ob- 
tained with the help of the plane-wave pseudopotential 
method for both the island facets and the wetting 
layer. Previous studies |l^,|l^ had missed the latter con- 
tribution, assuming that the surface energy of the wetting 
layer would be unchanged by the 3D transition. Our the- 
ory reproduces very well experimental data for the island 
size dependence on coverage. In particular, we improve 
over previous approaches by showing that the selectiv- 
ity of growth of a certain island size can be explained 
without invoking delicate elastic edge effects or island 
interactions . Furthermore, we demonstrate how the 
critical layer thickness depends on growth conditions, set- 
tling this long-standing issue. 

We propose a view of the growth process divided in 
three phases: an early nucleation phase which mainly de- 
termines the island density n, a second phase where the 
islands grow mostly on expense of the wetting layer, and 
a third phase characterized by Ostwald ripening. Since 
we are mostly interested in island sizes, we concentrate 




FIG. 1. Schematic illustration of the formation of coherent 
islands on the substrate surface. Oq and 9 are the nominal 
coverage and wetting layer thickness, respectively, a is the 
tilt angle of island facets and a is the island base length. 
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TABLE I. Surface energies 75 and surface stresses ctj,, ay 
for InAs surface reconstructions with the chemical potential 
MAs = MAs{buik) ~ 0.2eV. 
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on the second phase, and briefly discuss the third phase 
later. As long as the wetting layer acts as a source for ma- 
terial, existing nuclei will grow rapidly. Hereby the island 
density n remains constant 1^ . We treat it as an input to 
our model noting that the island density is determined by 
the growth kinetics. Furthermore, we assume the islands 
have identical shape and volume V. In the following, we 
discuss the island size in terms of a constrained thermo- 
dynamic equilibrium between islands and wetting layer, 
for a fixed island density. Although our approach is not 
limited to a certain system, here we consider, as an ex- 
ample, the strained {110} pyramidal shaped, dislocation- 
free InAs islands with a square base (area a^) on the 
GaAs(lOO) surface (with a wetting layer). We will also 
discuss the {111}/{111} faceted pyramidal islands later. 
We choose this system with a lattice mismatch Aa w 7%, 
because a large number of experimental and theoretical 
studies have been done j2|-|^,p. The real island shape 
may be more complex, but the simple island shape used 
here should still allow us to capture the important fea- 
tures of the island formation (see below and Ref. [p^). 
Figure 1 schematically illustrates the island formation on 
the substrate surface, and 9 are the nominal coverage 
and wetting layer thickness after island formation, re- 
spectively. We omit the interaction between islands (as 
we will discuss below this is a very good approximation in 
our case), as well the energies of edges and corners. The 
entropic contribution to free energy is also neglected. 

The total energy gain per unit volume of a single island 
can be expressed as 

Etot/V = et - ef,^ + [Sjf - a'^^^{B^)\|V 

+ (l/n - a") X [7wi((?) - 7wi(0o)]/^ (1) 

where eP} and e^jj^ are the elastic energy densities of is- 
land and uniformly strained film. The third term de- 
scribes the change in surface energy due to the island, 
with 7t being the surface energy of the island facets and 
S their area. The fourth term accounts for the thinning 
of the part of the wetting layer which feeds the island. 
7wi(^o), 7wi(^) are the formation energy of the wetting 
layer as a function of its thickness 6*, measured relative to 
InAs bulk kept at the GaAs lattice constant. This allows 
us to introduce the elastic contribution to the formation 
energy as a uniquely defined separate term egjjjj, similar 
to earlier work |14[| . From mass conservation, the volume 
of an island, V , is given by = ^a?" tana = ^(^^o ^ 




-o- 



^■"0 1 2 3 4 
coverage 5 (ML) 

FIG. 2. Formation energy of the wetting layer as a func- 
tion of thickness defined by 7wi^ = — MGaAs(buik)-^Ga 
-MinXl^buik)^!" - MAs(A^As - iVca - iVin), where A is the sur- 
face area and A*'; (i = As, Ga, In) are the number of particles of 
the species i in the supercell. E*°* is the total energy of the su- 
percell. /iGaAa(buik) and A'inAs(buik) are the chemical potential 
of GaAs bulk and of strained InAs bulk with the theoretical 
equilibrium lattice constant of GaAs bulk. From various con- 
figurations with A^in = 86 In atoms per (2 x 4) surface unit 
cell, those with the lowest formation energy are presented. 

where a and L are the tilt angle of island facets and 
the monolayer (ML) thickness, respectively. Equation 
(1) holds true as long as the island contains a sufficiently 
high number of atoms {e.g. 5 000 atoms) because elas- 
ticity theory is applicable and the reconstructions on the 
facets are completed [0. 

In order to obtain accurate values for the surface en- 
ergies and intrinsic surface stresses, these are computed 
using slab models of the surfaces with the help of the 
pseudopotential plane-wave method Q. After optimiz- 
ing the atomic geometries using consistently calculated 
forces on the atoms, the total energies of the slabs are 
computed, and the formation energies of various surfaces 
are obtained by subtracting the calculated total energy of 
an appropriate amount of bulk material. We further take 
into account the surface stress contribution to the surface 
energy up to the linear term |l^ , and a term proportional 
to the chemical potential /iAs of the environment in case 
of non-stoichiometric surfaces. Since epitaxial growth is 
mostly performed under As-rich conditions, all surface 
energies are evaluated close to equilibrium with bulk ar- 
senic {i.e. has = A'As(buik) ~ 0.2eV). For each facet, we 
have selected the reconstruction with the lowest energy 
from several candidates [ p^ . For the wetting layer, we 
consider the f32{2 x 4) reconstruction which is usually 
found on GaAs(OOl) and InAs(OOl) surfaces under mod- 
erately As-rich conditions ||l^. The results are given in 
Table I and Fig. 2. 

The elastic energy is calculated within continuum elas- 
ticity theory using the experimental elastic moduli to de- 
scribe the elastic properties of both the island and the 
substrate. For the island plus a 240 A thick slab rep- 
resenting the substrate, a finite-element approach is ap- 
plied to solve the nonlinear elasticity problem ||lj,|l6| . We 
calculate the elastic energy for several island shapes (dif- 
ferent tilt angles a of island facets) with a fixed island 
volume. For each particular shape, we can extract efl 
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FIG. 3. Total energy gain by islanding and various en- 
ergy contributions (solid lines) for n — lO^^cm^^, 6*0 = 1.8 
ML. The dashed line is the total energy gain for 
n = lO'^'^cm"^, ^0 ~ 1-5 ML. The arrows mark the minima 
of the curves. 

from the observed linear scaling relation with the island 
volume. The elastic energy density e^jj^^ of the uniformly 
strained film can be obtained by extrapolating the results 
for the islands to a=0, and we find it is in very good 
agreement with the value from linear elasticity theory. 

In Fig. 3, we show the various energy contributions 
and the total energy gain per volume for n = 10^^ 
cm^^j^o = 1.8 ML. The elastic relaxation energy (the 
first and second term in Eq.(l)) is negative due to strain 
relief and scales linearly with the island volume. Surface 
energy (the third term in Eq.(l)) is a cost, and therefore 
its contribution is positive. The wetting layer energy con- 
tribution (the fourth term in Eq.(l)) is also positive and 
depends complexly on the island volume, island density 
and coverage. We also show the energy contribution of 
the edges in Fig. 3, which becomes negligible compared 
to the other contributions for a large island (estimated as 
in Ref. [1^). It is important that an energy minimum ex- 
ists in the total energy gain curve. This indicates an op- 
timum island size can be obtained under certain growth 
conditions. The minimum in Fig. 3 corresponds to an 
island with about 38 000 atoms, which compares rea- 
sonably well with typical experimental values (between 
20 000 D and 50 000 atoms per island |l|]). The quite 
uniform islands prior to ripening observed in Ref. px| ] 
lend further support to the existence of an optimum is- 
land size. However, the island size strongly depends on 
the island density. Fig. 4(a) shows that the equilibrium is- 
land volume y is a hyperbolic function of island density 
n. As the experimental observations [p^ have shown, 
the island radius varies exponentially with the growth 
temperature. This was attributed to a kinetically lim- 
ited process It is consistent with the theory pre- 
sented here, since the island density is known to depend 
strongly on the growth temperature [ p^|j20|] , and thus ki- 
netics controls the growth through controlling the island 
nucleation density. Fig. 4(b) shows that our theory can 
reproduce very well the increase of the island radius with 
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FIG. 4. The dependence of the optimum island size (vol- 
ume and half-base) on the island density and the nominal 
coverage. In (a), the total energy gains for the optimum 
island size in various island densities are also shown. The 
experimental values in (b) are taken from Ref. [21] (•) and 
estimated from Ref. [3] (o). 



the amount of deposited material observed in various ex- 
periments, by using suitable island densities as input. 
The good agreement between theory and experiment also 
justifies the neglect of repulsive interactions between is- 
lands in the present study. A careful check indicates that 
the distances between the islands are quite large (larger 
than 65 nm and 220 nm for the high [2.3x10^° cm~^] 
and low [1.9x10^ cm~^] island densities, respectively). 



The island density 2.3 x 10 cm~ used to fit the exper- 
imental results |2l] agrees well with the experimentally 
estimated island density 1.5 — 2 x lO^'' cni^^. 

The critical layer thickness is a very interesting issue 
and the reported values vary from 1.2 to 2 ML [pT|-p3[. 
Our present theory puts us in position to discuss the criti- 
cal layer thickness, because the total energy gain depends 
on the coverage (see Fig. 3). When we deposit less ma- 
terial, keeping n fixed, the energy minimum rises above 
zero, i.e. island formation is no longer favorable. We 
take the critical layer thickness as the coverage at which 
the minimum energy equals zero (the error ± 0.01 ML). 
Our results, in Fig. 5, indicate that the critical layer thick- 
ness varies from 1.20 to 1.79 ML when the island density 
varies from 10^ cm~^ to 3.5x10^^ cm~^. Our theoreti- 
cal prediction matches the experimentally observed range 
(1.2-2 ML). We note that an exact determination is diffi- 
cult, because in the experiment various other factors may 
influence the critical layer thickness, like more complex 
island shapes, details of the growth method and growth 
conditions (e.g. III/V ratio) |^,^, or a possible corre- 
lation between island density and coverage ||^ . 

We have also performed a corresponding analysis for 
the strained InAs{lll}/{lil} pyramidal islands. The 
results show similar behavior as for the InAs{110} is- 
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FIG. 5. Critical layer thickness as a function of the island 
densities. 



lands and indicate that the basic features of our model 
do not depend on shape assumptions. The calculations 
also indicate that a larger island tends to be a steeper 
one, i.e. a {lll}/{ni} faceted pyramidal island, due 
to enhanced elastic energy relief |l^j2^. Our study with 
strained {lll}/{Tli} pyramidal islands also shows that 
the critical layer thickness somewhat depends on the is- 
land shape (see Fig. 5). However, we can still predict 
trends, e.g. for high growth temperatures (having a small 
island nucleation density) the critical layer is thinner. 
Our study also indicates that in case of a thinner critical 
layer, the island embryo should be larger than that for a 
thicker one. This can be understood in terms of a larger 
energy barrier which must be overcome by the embryo 
when nucleating on a thinner wetting layer. 

Finally, we briefly comment on the ripening of the is- 
lands. When no more material is supplied by the wetting 
layer, the island density is no longer constant, because 
smaller islands will dissolve again. Allowing the island 
density to vary, we find that larger islands at a lower den- 
sity are energetically preferred (see Fig. 4(a)). Thus, our 
theory is in accord with the observed Ostwald ripening 
pO| . However, since noticeable changes in the island size 
and density resulting from ripening typically take many 
days, ripening is not important for device applications, 
where the islands are covered by a capping layer after a 
very short growth interruption, and was not observed in 
previous experiments performed on a shorter time scale. 

In conclusion, we presented a novel explanation for the 
selection of particular sizes of self-assembled coherent is- 
lands in highly mismatched heteroepitaxy. It is possible 
to select the island size by changing the growth condi- 
tions and the nominal coverage. Our theory reproduces 
very well the experimental trends observed in the island 
growth. We have also shown how the critical layer thick- 
ness depends on growth conditions and settled this long- 
standing issue. 
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